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We discuss the influence of a uniform current, j , on the magnetization dynamics of a ferromagnetic 
metal. We find that the magnon energy t{q) has a current-induced contribution proportional to q-J , 
where J is the spin-current, and predict that collective dynamics will be more strongly damped at 
finite j. We obtain similar results for models with and without local moment participation in the 
magnetic order. For transition metal ferromagnets, we estimate that the uniform magnetic state 
will be destabilized for j > lO^A cm~^. We discuss the relationship of this effect to the spin-torque 
effects that alter magnetization dynamics in inhomogeneous magnetic systems. 

PACS numbers: 



INTRODUCTION 

The strong and robust magnetotransport effects that 
occur in metallic ferromagnets (anisotropic, tunnel, and 
giant magnetoresistance for example U) result from the 
sensitivity of magnetization orientation to external fields, 
combined with the strong magnetization-orientation de- 
pendent potentials felt by the current-carrying quasipar- 
ticles. This fundamentally interesting class of effects 
has been exploited in information storage technology for 
some time, and new variations continue to be discov- 
ered and explored . Attention has turned more recently 
to a distinct class of phenomena in which the relation- 
ship between quasiparticle and collective properties is in- 
verted, effects in which control of the quasiparticle state 
is used to manipulate collective properties rather than 
vice-versa. Of particular importance is the theoretical 
prediction of current induced magnetization switch- 
ing and related spin transfer effects in ferromagnetic mul- 
tilayers. The conditions necessary to achieve observable 
effects have been experimentally realized and the predic- 
tions of theory largely confirmed by a number of groups 
0,0,0,0,11,01 over the past several years. 

Current-induced switching is expected 0, 0| to oc- 
cur in magnetically inhomogeneous systems containing 
two or more weakly coupled magnetic layers. The work 
presented in the present paper was motivated by a re- 
lated theoretical prediction of Bazaily, Jones, and Zhang 
(hereafter BJZ), who argued that the energy functional of 
a uniform bulk half-metallic ferromagnet contains a term 
linear in the current of the quasiparticles [lH . i.e. that 
collective magnetic properties can be influenced by cur- 
rents even in a homogeneous bulk ferromagnetic metal. 
The current-induced term in the energy functional iden- 
tified by BJZ implies an additional contribution to the 
Landau-Lifshitz equations of motion and, in a quantum 
theory, to a change proportional to q ■ j in the magnon 



energy e{q). (Here qis the magnon or spin- wave wavevec- 
tor and j is the current density in the ferromagnet.) The 
BJZ theory predicts that a sufficiently large current den- 
sity will appreciably soften spin waves at finite wavevec- 
tors and eventually lead to an instability of a uniform 
ferromagnet. The current densities necessary to produce 
an instability were estimated by BJZ to be of order 10® A 
cm~^, roughly the same scale as the current densities at 
which spin-transfer phenomena are realized 00000, 
apparently suggesting to some that these two phenomena 
are deeply related. 

In this paper we establish that modification of spin- 
wave dynamics by current is a generic feature of all uni- 
form bulk metallic ferromagnets, not restricted to the 
half-metallic case considered by BJZ. We find that, in 
the general case, the extra term in the spin wave spec- 
trum 

Se{q) cx q- J (1) 

where is the spin current, i.e., the current carried by 
the majority carriers minus the current carried by the 
minority carriers the half metallic case J = j, 

recovering the result of Reference llj. For reasons that 
will become apparent later, we refer to the extra term in 
the spin wave dispersion as the spin wave Doppler shift, 
although this terminology ignores the role of underly- 
ing lattice as we shall explain. We also study the effect 
of a uniform current on spin- wave damping. The usual 
Gilbert damping law 7 cx e{q = 0), has an additional 
contribution proportional to the spin-current density. In 
our picture, a uniform current modifies collective mag- 
netization dynamics because it alters the distribution of 
quasiparticles in momentum space. 

Our paper is organized as follows. In Section II we 
present two general qualitative arguments which partially 
justify Eq.(^, independent of any detailed microscopic 
model. In Section III we substantiate the arguments with 
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a microscopic calculation of the spin wave spectrum for a 
ferromagnetic (but not necessarily half-metallic) phase of 
a Hubbard model, including the effect of the current. We 
derive Eq.(^, and demonstrate explicitly that when gen- 
eralized from the half-metallic case to the general case, 
the spin wave Doppler shift is proportional to the spin- 
current not the total current. The microscopic calcu- 
lation of Section III uses an effective action approach, 
which separates collective and quasiparticle coordinates 
in a natural way and is well suited to study their inter- 
play. In Section IV we specialize to the half-metallic case 
and re-derive the results of reference for the case of 
an s—d model ferromagnet. This serves the purpose of es- 
tablishing a clear formal connection between the deriva- 
tions presented in Sections II and III and the derivation 
presented by BJZ, which appear superficially to be quite 
distinct. In Section V we discuss the effect of a current 
on spin wave damping. We consider both damping due 
to the coupling of spin waves with the quasiparticles and 
two magnon damping, which we argue is enhanced by 
the spin wave Doppler shift of Eq.(^). In Section VI we 
discuss the relationship between the spin wave Doppler 
shift and spin-transfer in inhomogeneous ferromagnets. 
Finally, in Section VII we summarize our main results 
and present our conclusions. 



QUALITATIVE EXPLANATION OF THE 
CURRENT-INDUCED MAGNON ENERGY 
SHIFT 

The low energy collective dynamics of the magneti- 
zation orientation in a ferromagnet is described by the 
Landau-Lifshitz equation: 



hdn{r,t) 
dt 



= f7 X 



sn 



_dn{f,t) 
* dt 



(2) 



where f2(r, t) is an unimodular vector field which de- 
scribes the orientation of the collective magnetization 

and E ^f2, 9if2j^ is an energy functional of <) and 
its derivatives. The generic applicability of this equa- 
tion follows from the collective nature of spin-dynamics 
in ferromagnets. It can be derived from a number of 
different microscopic models in a number of different 
ways. In particular, this equation describes the low- 
energy long-wavelength dynamics of the two models of 
metallic ferromagnetism that we consider in later sec- 
tions. Normally E is minimized by a coUinear configura- 
tions fi(r, t) = f^o along some privileged easy direction. 
The Landau-Lifshitz equations linearized around r2o have 
solutions which describe distortions of the magnetization 
orientation that propagate like waves with wave q and 
frequency tjj{q) [l3|- In a quantum treatment, magneti- 
zation orientation fluctuations are quantized in units of 



In a metallic ferromagnet, the quasiparticles occupy 
bands [13 that are energetically split by an effective 
Zeeman-coupling magnetic field oriented along the direc- 
tion Non-coUinear configurations are penalized be- 
cause band-electron kinetic energies are raised by an in- 
homogeneous effective field Q{f,t). The easy axis is de- 
termined by spin-orbit interactions of the band electrons 
and by the magnetostatic energy, which because of its 
long range depends on the overall shape of the sample. 

The dynamics generated by the first term in square 
brackets in Eq. is energy conserving whereas the sec- 
ond term, proportional to the dimensionless coefficient a, 
transfers energy from the collective coordinate to other 
degrees of freedom. In a metallic ferromagnet, the damp- 
ing is partly due to the excitation of electron-hole pairs 
in response to the temporal evolution of Q. It is clear, 
therefore, that there is an intimate relation between the 
dynamics of the collective coordinate and the state of the 
quasiparticles. What's more, when current flows inside a 
ferromagnet, the momentum-space distribution functions 
that describe quasiparticle state occupation probabilities 
are altered. It is natural, therefore, to expect that the 
dissipative dynamics of the collective magnetization will 
be affected by current flow. In Ref. 0it was shown that, 
in a half metallic ferromagnet modeled by a s-d model (a 
model with a single band coupled by exchange interac- 
tions to local moments), the energy functional E has a 
term linearly proportional to the quasiparticle current, j. 
In the following paragraphs we present three arguments 
to support the idea that the spin wave spectrum of any 
metallic ferromagnet is modified by a uniform current in 
a manner similar to that suggested by Eq. ( ^ . 

We start with the simplest half-metallic ferro- 

magnetic electron gas, in which the current effect can be 
understood simply in terms of Galilean invariance. The 
current carrying state of this system is simply one in 
which the entire electronic systems moves along with a 
drift velocity vd- A spin wave excitation is one in which 
the magnetization orientation precesses around the easy 
axis with frequency uj{q): 

fl = (e sin((f- r — wo(9)t), e cos(q- f— ujQ{q)t), 1 — e/2) (3) 

In the lab frame, the system is seen as moving with ve- 
locity t//), and carrying current j = —nevo- The fixed 
position r^ in the lab frame, has position "Tl — vot in 
the moving frame. The precession frequency seen at a 
fixed lab frame position is therefore Doppler shifted to 
^0(9) + 9 ■ VD- 

This simple effect is the essence of the spin-wave 
Doppler shift. In terms of the current density the spin- 
wave Doppler shift in the magnon energy is hq ■ j /en. 
Systems of practical interest are neither Galilean invari- 
ant nor, with a few possible exceptions, half metallic, 
however so a more detailed analysis is required to deter- 
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mine how the spin-wave Doppler shift is manifested in 
real systems. 

A second useful point of view follows from considering 
a single-mode-approximation for the quantum spin- wave 
energy e(g) = hu!{q). Elementary magnon excitations of 
a ferromagnet reduce the total spin projection along the 
easy axis by one unit and add crystal momentum hq. A 
state with the correct quantum numbers can be generated 
starting from the ferromagnetic ground state (or from a 
state that carries a uniform current) |^'o) by acting on it 
with the 'magnon creation operator' 



E 

i=l,N 



s_i exp(i(f ■ 



(4) 



where S-i is the spin-lowering operator for the z-th parti- 
cle. Two-particle Greens functions constructed from this 
operator have poles with large residues at magnon excita- 
tion energies. The single- mode approximation consists of 
using \'^{q)) = s^{—q)\'i'o) as a variational wavefunction 
for the magnon state at wavevector q. Given this ap- 
proximation for the magnon state, its excitation energy 



parabolic bands, in generalizing the magnon Doppler ef- 
fect from half-metallic ferromagnets to ferromagnets with 
states of both spins occupied, the current is replaced by 
the spin-current ^7, and the density by the spin-density. 

Finally, the same result can be derived by considering 
a variational wave function for the spin-wave state of a 
ferromagnetic metal in which all quasiparticle states that 
are singly occupied share a common spinor that describes 
long- wavelength spatial precession around the easy direc- 
tion. For example if the x direction is the easy direction 
the spinor that describes small amplitude precession is 
{u,v) = (1 — 77^/2, ryexp(ig • r)). The q ■ J correction 
then follows by observing that the magnon energy equals 
the energy change divided by the change in the x direc- 
tion magnetization component, with both quantities be- 
ing proportional to rf' at small r\. This findings suggest 
that the explicit approximate expression for the magnon 
Doppler shift, derived from the SMA for parabolic bands, 
is likely to qualitatively correct even for realistic ferro- 
magnets with more complicated band structures. Indeed, 
that is the conclusion that follows from the more micro- 
scopic derivations in the following two sections. 



(5) 



can be expressed in terms of the expectation value of a 
commutator between the general many-particle Hamilto- 
nian H and either magnon creation or annihilation oper- 
ators and simplified to the following form: 



hq 
m 



h^q^ 
— — + 
2m 

Ej,(*o|s+iS_jexp[jg- {t 



(«'o|s+(q)s-(-(7)|1'o) 



(6) 



The second term on the right hand side of Eq.( O is 
the magnon Doppler shift term. In this term s±i and pi 
are the spin raising and lowering and momentum oper- 
ators for particle i. The numerator and denominator of 
this term are, in general, complex two-particle correla- 
tion functions. The correlation functions are simplified 
when the ferromagnetic state is approximated by a Slater 
determinant with definite occupation numbers for both 
majority {]) and minority (J,) spin momentum states, i.e. 
by the electron gas Stoner model ferromagnetic ground 
state. Then to leading order in q we find that the magnon 
Doppler shift has the value 



Seiq)^"^-- 

m n| 



J 



(7) 



Eq. Ul is most easily obtained by writing the operators 
whose expectation values need to be evaluated as a sum 
of one-body and two-body terms and then using standard 
second quantization identities. The most important con- 
clusion suggested by this equation is that, at least for 



CURRENT DRIVEN SPIN WAVES IN A 
HUBBARD MODEL FERROMAGNET 

In order to explain our theory of the influence of uni- 
form currents on the spin- wave spectrum, we first recall 
how spin-waves and quasiparticle states are related in 
equilibrium. This development will also establish the no- 
tation we use for the non-equilibrium case. The descrip- 
tion we use is one in which a collective fluctuation field 
interacts with fermionic quasiparticle fields. It allows us 
to borrow from standard theories of quantum harmonic 
oscillators weakly coupled to a bath, in order to general- 
ize the theory of collective dynamics from equilibrium to 
non-equilibrium cases. 



Hamiltonian and effective action 

In the previous section we discussed three general argu- 
ments in support of the existence of a spin wave Doppler 
shift in a metallic ferromagnet that is proportional to the 
spin current as in Eq.(^. We now look more closely at 
the underlying physics by carrying out an explicit mi- 
croscopic calculation of the spin waves for a Hubbard 
model in the presence of a current. Unlike the s-d model 
considered in Ref. [ill the Hubbard model allows for fer- 
romagnetism in a system with only itinerant electrons. 
The Hubbard model Hamiltonian is 0| : 



(8) 



The elementary excitations of a metallic ferromagnet 
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are quasiparticles and spin waves. We want to derive the 
propagator for the spin waves of the ferromagnetic phase 
of this model and to see how is affected by a quasiparticfe 
current. To do so, it is convenient to use the functional 
integration approach, 0,0,0,0], in which the quasi- 
particles are integrated out and an effective action for the 
spin waves is obtained. This procedure is sketched be- 
low, the details can be found in Refs. 
The final result for the spin wave spectrum is equivalent 
to that obtained by doing a Random Phase Approxima- 
tion (RPA) calculation. However, the effective action 
approach provides a convenient conceptual framework to 
understand the connection between spin waves and non- 
equilibrium quasiparticle states, the central focus of this 
paper. 

The interaction term in the Hubbard model can be 
written as 



u 



We represent the partition function of this model as a 
path integral over fcrmion coherent states labeled 
by {^Q, , where a = i,a. The key idea which allows 
quasiparticle and collective degrees of freedom to be sep- 
arated, while still treating the magnetization as a quan- 
tum field, is the introduction of a Hubbard-Stratonovich 
transformation (ilj to represent the interaction term. By 
making this transformation, we trade a problem of inter- 
acting fermions for a problem of independent fermions 
whose spin is coupled to a bosonic spin-splitting effective 
magnetic field Ai(r), which acts as the collective mag- 
netic coordinate. The partition fimction reads: 



I?*„(r)I?^'„(r)I?A,(T)exp 

where the action is 

_ 3A,(r)2 



Jo 



m 



Y.^^^-i^)Sr,la''^jAr) (9) 



and 



u 



-5.,, (10) 



is the inverse of the Green's function operator. 

The action @ is the sum of three terms, i) non- 
interacting tight binding fermions (with a Hartrcc shift), 
a) a term quadratic in the bosonic field and in) their 

coupling • Si, where Si[t) = J2a,cr' l'^i,crTa,cr'^i',a'- 

Since the action is quadratic in the fermion variables, 
the fermion functional integral can be formally evalu- 
ated. This allows to write the partition function as a 
path integral over the auxiliary field A only, 



Z = / VA,iT)e 



-Sett (A) 



(11) 



where the effective action reads 

'■f' 3A,(r)^ 
8U 



5eff(A)= / drV 



Tr Ln 



g-\A) (12) 



Eas.flll|) and f ll2|) are one of the many possible represen- 
tations of the exact partition function for the Hubbard 
Model. The effective action (I12|l describes a complicated 
quantum field theory for Ai(r). 



Mean Field theory: Spin-split bands 

The first step in a field theory of ferromagnetism is 
usually to look for classical solutions, i.e. for field config- 
uration Af (r) for which the effective action is stationary. 
The saddle point equation reads Af = ^{Si), where 
the average is computed with a Green function ^(A'^') 
obtained by replacing, in Eq.(^J), the fluctuating field 
Ai (r) by the saddle point solution. 

Assuming the existence of a ferromagnetic mean-field 
state, the classical solution for a perfect crystal is static 
(independent of r) and homogeneous (independent of i). 
It is therefore characterized by a direction n and a length 
jA"^'! = A. Because of the spin rotational invariance of 
the Hubbard Hamiltonian, n is arbitrary. In real systems 
n is determined by spin-orbit interactions and magneto- 
static effects. The mean-field Green's function, C/(A'^'), 
describes fermions which occupy bands that are spin-split 
by an effective magnetic field along n (See Fig.(l)). The 
magnitude of the spin splitting. A, is obtained from the 
saddle point equations, which, for this simple model, re- 
duce to the following form: 



^ 4C/ 1 r 



(13) 



where — e(fc) — are the quasiparticle energies of 
the spin-split bands and M is the number of lattice sites. 
Notice that the majority band has spins parallel to n, de- 
noted by t- The saddle point equations show explicitly 
that the auxiliary field A'^' is proportional to the aver- 
age fermion magnetization, which usually appears as the 
fundamental field in classical micromagnetic theories for 
realistic magnetic materials. Hereafter we refer to A(t) 
as the 



Spin waves without current 

We are interested in the dynamics of the collective co- 
ordinate, so that the static solution obtained by solving 
the mean field approximation is insufficient. To describe 
the elementary collective excitations we need to study 
small amplitude dynamic fluctuations of the collective 
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k (a.u.) 



FIG. 1: Mean field quasiparticle bands. Dashed line shows 
the Fermi Energy. A is the spin splitting energy. 



coordinate around the static solution: 



A.(r) 



(14) 



We introduce Eq.(^^ into the effective action (Eq.(El)) 

r ^ -1 3 

and neglect terms of order SAi{T) and higher. The 

resulting actionS'ci(A'^') + 5sw, where the first term is 
the classical approximation to the effective action and 
the fluctuation correction is: 



'S'sw 



1 



Q,a,b 



V SAa{Q)lCab{Q)SAt{-Q) (15) 



where Q is a shorthand for q, ivn, and a, b stand for Carte- 
sian coordinates. Note that the bosonic fields, SA{Q) are 
dimensionless and the Kernel /C has dimensions of inverse 
energy. This action defines a field theory for the spin 
fluctuations. The equilibrium Matsubara Green function, 
Vabiq, ifn) , is given [i^,!!^ by the inverse of spin fluctu- 
ation Kernel, JCabiQ)- Analytical expressions for ICab{Q) 
are readily evaluated for the case of parabolic bands and 
are appealed to below. We obtain the retarded spin fluc- 
tuation propagator by analytical continuation of the Mat- 
subara propagator: D^°i^{q,uj) = Vab^q^ivn — > w + iO+) 
The imaginary part of the retarded propagator summa- 
rizes the spectrum and the damping of the spin fluctua- 
tions most directly. 

The theory defined by Eq. ( I15|l includes two types of 
spin fluctuations which are very different: i) longitudi- 
nal fluctuations (parallel to n), or amplitude modes and 
ii) transverse fluctuations (perpendicular to n), or spin 
waves. The amplitude modes involve a change in the 
magnitude of the local spin splitting. A, and are either 
over damped or appear at energies above the continuum 
of spin-diagonal particle-hole excitations. In contrast, 
the spin waves are gapless in the limit g = 0, in agree- 
ment with the Goldstone theorem, and are often weakly 
damped even in realistic situations. Note that the ampli- 
tude modes decouple from the spin wave modes for small 



amplitude fluctuations. For x 

ICabiQ) = 





n, we can write 







-yy 



y^ 



(16) 



Since the low energy dynamics of a metallic ferromag- 
net is governed by transverse spin fluctuations, we do 
not discuss longitudinal fluctuations further. After ana- 
lytic continuation, we obtain the following result for the 
inverse of the retarded transverse spin fluctuation Green 
function (_D''°')~^, which is diagonal when we rotate from 
y, z to +z ± iy chiral representations. The diagonal ele- 
ments are then 

4U 1 



3 1 + |[/F(±g,±w) 



(17) 



where T(q,uj) is the Lindhard function evaluated with 
the spin-split mean-field bands: 



ri- 
fe 



k+q 



k k 



k+q 



^0^ 



(18) 



where n'l is shorthand for the Fermi-Dirac occupation 
function np 
bers. Eqs. 



for the quasiparticle occupation num- 

( I17|l and ( I18|l make it clear that the spin 
wave spectrum is a functional of the occupation function 
np for the quasi-particlcs in the spin-split bands. The 
infiuence of a current on the spin- wave spectrum will en- 
ter our theory through non-equilibrium values of these 
occupation numbers. 

In the case of parabolic bands (still without current), 
the Taylor expansion of the Lindhardt function in the 
low-energy low-frequency limit gives the following result 
for the spin wave propagator: 

6 uj ± pq'^ 

where p is the spin stiffness which is easily computed an- 
alytically for the case of parabolic bands. The poles of 
Ea.f ll9|l give the well known result for the spin wave dis- 
persion, UJ ~ ipq^. Several remarks are in order: i) In 
real systems, spin-orbit interactions lift spin rotational 
invariance, resulting in a gap for the q = spin waves. 
The size of the gap is typically of order of 1 /leV (24j . 
ii) The interplay between disorder and spin orbit inter- 
actions, absent in the above model, gives rise to a broad- 
ening of the spin wave spectrum, even at small frequency 
and momentum. In Section V we address this issue and 
discuss how damping is changed in the presence of a cur- 
rent. 



Spin waves with current 

In the previous subsection we derived the spin wave 
spectrum of a metallic ferromagnet in thermal equilib- 
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rium. Equations 1)17(1 and (|18() establish a clear con- 
nection between spin waves and quasiparticle distribu- 
tions. In order to address the same problem in the pres- 
ence of a current, a non-equilibrium formalism is needed. 
By taking advantage of the formulation discussed above 
in which collective excitations interact with fermion 
particle-hole excitations we are able to appeal to estab- 
lished results for harmonic oscillators weakly coupled to 
a bath. In the equilibrium case, the fact that the low- 
energy Hamiltonian for magnetization-orientation fluctu- 
ations is that of a harmonic oscillator follows by expand- 
ing the fluctuation action to leading order in lo to show 
that y and z direction fluctuations are canonically con- 
jugate. In our model magnons are coupled to a bath 
of spin-flip particle-hole excitations. Following system- 
bath weak coupling master equation analyses '23 we find 
that the collective dynamics in the presence of a non- 
equilibrium current-carrying quasiparticle system differs 
from the equilibrium one simply by replacing Fermi occu- 
pation numbers by the non-equilibrium occupation num- 
bers of the current-carrying state. The following term 
therefore appears in the Taylor expansion of the Lind- 
hardt function F: 



dV 



1 



g— cj— 



AAA- 



E 



de{k) 



T 
k 



k 



(20) 



Since this expression uses the easy direction x as the spin- 
quantization axis, the x (spin) component of the spin 
current is: 



e de{k) 



dk 



T i 
ni — nt 
k k 



(21) 



case, however, the explicit expression was derived for 
the case of free-particle parabolic bands only. Eq. ( 
states that the spin wave spectrum of metallic ferromag- 
net driven by a current is modifled in proportion to the 
resulting spin current. 

In the half metallic case, when the density of minority 
carriers is zero, the spin current is equal to the total 
current and we recover the result of BJZ flll|. In that 
limit A = and p ~ leading to 



Lu ^ —q q-j^—q -hq-Vo 

2m en Zm 



(24) 



where we have expressed the current as j — envo with 
vd the dri/t velocity, generalizing the half-metallic simple 
Doppler shift result to non-parabolic bands. 



Spin wave instability 

Eqs. { \2'6\ and f I24() . taken at face value, predict that 
the energy of a spin waves is negative and therefore that 
the uniform ferromagnetic state is destabilized by an ar- 
bitrarily small current. If this were really true, it would 
presumably be a rather obvious and well known experi- 
mental fact. It is not true because spin waves in real fer- 
romagnetic materials have a gap due to both spin-orbit 
interactions and magnetostatic energy. Inserting by hand 
this (ferromagnetic resonance) gap, the spin wave disper- 
sion reads: 



UJ — ujo + pq 



2U h_ 
3A 



-q-J 



(25) 



so that 



dr 

dqi 



(22) 



The quantity J7i, the component of the spin current po- 
larized along the magnetization direction w = x and 
flowing along the i axis, is the difference between the 
current carried by majority and minority carriers. In 
equilibrium there is no current and no linear term oc- 
curs in the wavevector Taylor series expansion, leading 
to quadratic magnon dispersion as obtained in Ea.f I19|l . 
When (charge) current flows through the ferromagnet, 
the difference in carrier density and mobility between 
majority and minority bands inevitably gives rise to a 
nonzero spin current 30] . We therefore obtain the follow- 
ing spectrum for spin waves in the presence of a current: 



P9 



3A e 



(23) 



This equation is the central result of our paper. No- 
tice that it is in precise agreement with the single-mode- 
approximation expression since A = "^{n^ —n\)\ in that 



so that it takes a critical spin current to close the spin 
wave gap. In Fig. (2) we plot the current driven spin wave 
spectrum assuming ujq = IfieV, the electronic density of 
iron (n = 1.17 10^'^ cm^^) and a Doppler shift given by 
q Vd. The critical current so estimated is ~ 1.1 10^ A 
cm~^ for a typical system. This critical current could be 
much lower, perhaps by several orders of magnitude, in 
metallic ferromagnets in which material parameters have 
been tuned to minimize the spin- wave gap. Experimental 
searches for current-driven anomalies in permalloy thin 
films, for example, could prove to be fruitful. 



Spin wave action with current 

In the small lu and small q limit, the spin waves are 
independent and their action is equivalent to that of an 
ensemble of non interacting harmonic oscillators, indexed 
with the label q. The Matsubara action for a single os- 
cillator mode is the frequency sum of 



[Pq J 




(26) 
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FIG. 2: Current modified spin-wave spectrum 



where the diagonal terms are the Hamihonian part of 
the action and the off-diagonal term can be interpreted 
as a Berry phase. For the spin waves, the analog of p 
and X are, modulo some constants, 6Ay,6Az. In this 
representation, the low to and low q spin wave action 
reads: 



pq ■ q —ILO 
iu) pq - q 



2U h ^ 



-i 

1 



(27) 



This representation makes it clear that the spin wave 
Doppler shift appears as a modification of the term 
which couples the canonically conjugate variables, SAy 
and SAz, i.e., the spin wave Doppler shift modifies the 
Berry phase. When expressed in this way, the spin-wave 
Doppler shift is partly analogous to the change in su- 
perfluid velocity in a superfluid that carries a finite mass 
current, and the stability limit we have discussed is partly 
analogous to the Landau criterion for the critical veloc- 
ity of a superfluid. These analogies are closer in the case 
of ideal easy-plane ferromagnets, which like superfluids 
have collective modes with linear dispersion instead of a 
having a gap. 



ALTERNATE DERIVATION OF SPIN- WAVE 
DOPPLER SHIFT 

In the previous section we have used a functional in- 
tegral approach to calculate how the spin wave propaga- 
tor of a Hubbard model metallic ferromagnet is modified 
when current flows through the system. The BJZ deriva- 
tion of the same effect was based on an identity at the 
operator level. BJZ used a, s — d model, i.e. a Hamilto- 
nian for itinerant (s) electrons interacting with localized 
(d) spins via an exchange interaction. They considered 
the limit of very large exchange interaction and low s 
electron density, so that, in the ground state, the elec- 
trons are fully spin polarized. They then introduced a 
local spin rotation transformation defined so that at ev- 
ery point of the space the spin of the s electrons is parallel 
with the local value of the d electrons magnetic moment. 
This unitary transformation has been previously used for 



both s — d and other microscopic models of ferromag- 
netism [H HI HHH [111 . In the transformed frame, the 
exchange interaction is always diagonal in the spin index 
but the expression for the kinetic energy is complicated, 
and includes new terms. One of the new terms couples 
the s electron current to a space derivative of the local 
spin magnetization. It is from this term in the exchange 
energy that BJZ derived the modification of the Landau- 
Lifshitz equations that we have identified as a spin- wave 
Doppler shift. In this section we bridge the gap between 
the two derivations. We recover the half metallic s-d 
Hamiltonian result of BJZ in a systematic way. 

The continuum s-d model describes itinerant electrons, 
"00-, interacting with a continuum of localized quantum 
spins, M{x), through a exchange interaction of strength 
J. The Hamiltonian for parabolic bands is given by: 



n 



2m 



M{x) 



where r are the Pauli matrices. In order to derive an 
effective theory for the collective behavior of this system, 
we express its partition function as a coherent state path 
integral: 



Z = 



-H 



where r is imaginary time, VL{x,t) — ^M{x,t) = 
[cos((/)) sin(6'), sin((/)) sin(6'), cos(6')] is the unimodular vec- 
tor field which labels the spin coherent states, Sb is the 
Berry phase term that captures the spin commutation 
relations '23'|, and ^' are the Grassmann numbers which 
label the fcrmion coherent states 22] . 

Following BJZ, we perform a unitary transformation 
on the spins of the itinerant electrons so that, at each 
point of time and space, the quantization axis is paral- 
lel to AI{x,t). BJZ considered only the limit of very 
strong ferromagnetic J, so that the spins of the occupied 
electronic states are always parallel to M{x,t) and we 
can ignore the antiparallel electrons. This approxima- 
tion is valid in half-metallic systems for energies much 
smaller than J, the local spin splitting. In this approx- 
imation the action for the parallel fermions in the ro- 
tated frame, denoted by $(a?, r), can then be written as 
S — Sb + S'o + -I- 5*2 where: 



Sb — 

50 = 

51 = 



f ^ f \ 

dr / dx i S cs { 1 ^ — cos(0)9r(/' 

J V 25cs/ 

dr dx^idr-p- -.JS\^ 



S2 — — / dr I dx 



Jp + Jd 



Am, (28) 
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cs is the density of local moments with spin S, and 



A{n) 



he 



cos(0) is an effective vector potential which 



depends on the local spin configuration, In Eq.EHl Jp 
and J7b are respectively paramagnetic and diamagnetic 
contributions to the current density defined by 

Jp = f¥(f,r)V$(f,T) - V$(f,T)$(f,^ 

e 



Jd 



(29) 



The Jp ■ A coupling has the form anticipated by BJZ. 
To address the magnetic elementary excitation spectrum 
we formally integrate out the fermion fields $ and expand 
to quadratic order in magnetic fluctuations. The action 
expressed in terms of only the spin fields is Scs{^) — 
5b + Tr[ln5~i] with 



g-\e,4>)^dr- 

+i cos{6)dr(l) H 



2m \ 2c 



(30) 



Expanding around the x {0 — ^ (p = 0) direction we 
obtain for the spin-wave action 



TrLn [g-\e,(b)] = TrLn \g-\-,o) + dg-\n) 

To leading order in 60 the action reads: 



ec 



dx 



j + — A{n) -Ain) 

mc 



(31) 



where j = Tr 



is the average current 



g(0,f)Jp(f,T) 

and n = Tr [G{0, ^)^{x,t)^{x,t)\ is the average den- 
sity in the coUinear ground state. In deriving this ex- 
pression we allowed the mean-field fermion quasiparticle 
occupation numbers to assume values consistent with a 
non-equilibrium current-carrying state. 

Equation defines a theory for the collective mag- 
netization of the ferromagnet. The first two terms are 
the Berry phase of the d spin and a renormalization of 
the Berry phase due to the spin of the s electrons, similar 
to that derived by Millis et al. for the double exchange 
model The third term describes the energy penalty 
for non coUinear configurations, or spin stiffness. The 
terms in the second line yield the coupling of the average 
(paramagnetic and diamagnetic) currents to the collec- 
tive magnetization. 

The semiclassical equations of motion of (13 111 yield the 
Landau Lifshitz (LL) equations including the jiViil x 
term derived by BJZ (equation (5)). In the case of BJZ, 



the LL equations are derived from a micromagnetic en- 
ergy functional plus the paramagnetic current term. In 
our case, the whole functional is derived from the micro- 
scopic Hamiltonian. The spin wave expansion for H31|l 
around a classical homogeneous ground state, fici = x 
is obtained by expanding = fici + <^^^ and identifying 
6 Vly ^ (j), 5 Q,z — cos{6). Dropping terms of order 5Vl^ 
and higher, the action (|3HI becomes: 



213V 



(32) 



Q,ab 



as in Eq. ( I15|l . After analytical continuation, the spin 
wave kernel 1, in the y, z representation: 



Cs 



pq 



-iS'uj 



• CI/ 2 

lb UJ pq 



+ -3 -q 

e 



~i 

1 



(33) 



where r = p = r j^, and S' = S + ^. The main 

cs ' " 4:m ' 2 

difference between s — d and Hubbard model result is 
the appearance here of both local moment and itinerant 
electron (r/2) contributions to the Berry phase, which is 
proportional to the total spin density. Note that since 
A is quadratic in the spin wave variables, the term 
in gives no contribution to (|32l) . After diagonaliza- 
tion of Eq. H33I) we obtain the retarded propagator for 
the spin wave variables. The real and imaginary part of 
the poles of the retarded propagator give the spin wave 
dispersion and damping, respectively. In this theory, the 
imaginary part is zero, since the spin flip of quasiparticles 
is blocked. The real part reads: 



h^ruf' 
'ASr 



2S, 



-J ■ q 



C5 1 



2Scs 



(34) 



Hence, we see how the spin wave dispersion in this 
theory has the q ■ j term derived by BJZ. Since the sys- 
tem described by the theory is fully polarized, the cur- 
rent and the spin current (polarized along the ground 
state magnetization direction) are identical. This result 
is to be compared with Eq. H24() , derived with a different 
method for a different microscopic model. We conclude 
that spin-wave Doppler shifts due to spin currents are 
generic, although their quantitative details can depend 
on the microscopic physics of the ferromagnet. 



ENHANCED SPIN- WAVE DAMPING AT FINITE 
CURRENT 

In Sections III and IV we have shown how the disper- 
sion of spin waves in a metallic ferromagnet is affected 
by current flow, and we have obtained results compatible 
with those of BJZ ■ this section we address a prob- 
lem which, to our knowledge, has remained unexplored 
so far: how does the current flow affect the lifetime of the 
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spin waves. In subsection A we analyze the damping of 
spin waves at zero current. In the subsections B and C 
we discuss how these results are modified by the presence 
of a current. 

A ferromagnetic resonance (FMR) experiment probes 
the dynamics of the coherent or q — spin wave mode. 
The signal lincwidth is inversely proportional to the co- 
herent mode lifetime, the time that it takes for a trans- 
verse magnetic fluctuation to relax back to zero. Spin 
waves have a finite lifetime because they are coupled to 
each other and to other degrees of freedom, including 
phonons and electronic quasiparticles. In ferromagnetic 
metals, the quasiparticles are an important part of the 
dissipative environment of the spin waves [sH Is^ Is^ Is^ l . 
and we can therefore expect that quasiparticle current 
flow affects the spin wave lifetime to some degree. In or- 
der to discuss this effect, it is useful to first develop the 
theory of quasiparticle spin- wave damping in equilibrium. 



Damping at zero current 

The elementary excitation energies for the ferromag- 
netic phase of the Hubbard model, are specified by the 
locations of poles in Ea.f I17|l . The damping rate is pro- 
portional to the imaginary part of the transverse fluc- 
tuation propagator. According to Ea.f I17|l . the damp- 
ing of a spin wave with frequency lu and momentum q, 
= -2Im [r(w, q)] is given by: 



7((f,w) 



27r 

k 



[4 



fc4 



4 



4 

fc4 



(35) 



In the absence of disorder, this quantity is nonzero when 
\q\ is comparable to kp^ — kpi or when w ~ A, the band 
spin-splitting. Either disorder, which breaks transla- 
tional symmetry leading to violations of momentum con- 
servation selection rules, or spin-orbit interactions, which 
cause all quasiparticles to have mixed spin character, will 
lead to a finite electronic damping rate at characteristic 
collective motion frequencies. Because this damping is 
extrinsic, however, its numerical value is usually diffi- 
cult to estimate. It is often not known whether coupling 
to electronic quasiparticles, phonons, or other degrees of 
freedom dominates the damping. 

Formally generalizing Eq. ( I35|l to the case with disor- 
der and spin orbit interactions leads to 



7(w) cx ^ S^-, 

k.k' ,u.u' 

(36) 

where Su^u>{k,k') = \{k,v\S^~^\k' ,v')\'^ is a matrix ele- 
ment between disorder broadened initial and final quasi- 
particle states, labeled by momentum k and band index 
V (but not Bloch states). Averaging out the extrinsic de- 
pendence on wavevector labels by letting Svy{k,k') 



Si,y we obtain 

7(w) = n^Y.^-y J '^^ J de'N,{e)N,,{e')x 

X {nie)-n{e'))5[e-e' +uj] (37) 

where Ni,{e) is the density of states of the band ly. For 
uj of the order of the ferromagnetic resonance frequency, 
we can expand Eq. 1)3 7|) to lowest order in uj: 



7(w) 



n2^5,,,-A^,(ef)iV,-(ef) 



(38) 



This result can be considered a microscopic justifica- 
tion of the Gilbert damping law, which states that the 
damping rate is linearly proportional to the resonance 
frequency and vanishes at w — 0. The proportionality 
between frequency and damping rate follows from phase 
space considerations: the higher the spin wave frequency 
LU, the larger the number of quasiparticle spin flip pro- 
cesses compatible with energy conservation. 



Damping at finite current 

We analyze how a current modifies quasiparticle damp- 
ing, we again appeal to the picture of magnons as har- 
monic oscillators coupled to a bath of particle-hole exci- 
tations and borrow results from master equation results 
for oscillators weakly coupled to a bath 25'| For magne- 
tization in the 'f direction, magnon creation is accom- 
panied by quasiparticle-spin raising and magnon annihi- 
lation is accompanied by quasiparticle-spin lowering. It 
turns out that only the difference between the rate 
of quasiparticle up-to-down and quasiparticle down-to-up 
transitions enters the equation that describes the magne- 
tization evolution. This transition rate difference leads to 
the same combination of quasiparticle occupation num- 
bers as in Eg. ( 1381) . except that the occupation numbers 
characterize the current-carrying state and are not Fermi 
factors. For metals we can use the standard approximate 
form|35l| for the quasiparticle distribution function in a 
current carrying state: 



dn 




^ 97 


fc-l 



(39) 



Because of the independent sums over k and k' in 
Eg.f I36|l . and because it is a simple difference of Fermi 
factors that enters the damping expression, we conclude 
that the quasiparticle damping correction will vanish to 
leading order in the spin-dependent drift velocities vf^. 
We reach this conclusion even though the phase space 
for spin-flip quasiparticle transitions at the spin-wave en- 
ergy is altered by a factor 1 when ep x ^ ~ eoj where 
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Ci? is a characteristic quasiparticle energy scale, i.e. the 
up-to-down and down-to-up transition rates change sig- 
nificantly when this condition is met, but not their differ- 
ence. To obtain a crude estimate for the current at which 
this condition is satisfied we use the following data|35j for 
iron: n « 1.7 10^^, Fermi velocity ~ 1.98 10^ cms~^ The 
drift velocity corresponding to a current density of lO^^A 
cm~^ is Vd = — 10^"^ cm s^^. The typical energy 
of a long-wavelength magnon is ^ 10"^ eV. Therefore, 
current densities of the order of 10^ A cm~^ and larger 
will substantially change the coupling of spin-waves to 
their quasiparticle environment. Although this change 
will influence the spin-wave density-matrix, magnetiza- 
tion fluctuation damping itself will not be altered by this 
mechanism until much stronger currents are reached. 



Two magnon damping 

In the previous subsections we have calculated the 
damping of the lowest energy spin wave due to its cou- 
pling to the reservoir of quasiparticles. In this section 
we study damping of the coherent rotation mode [q — 
spin wave) due to its coupling to finite q spin waves. This 
mechanism is known as two magnon scattering and is ef- 
ficient when the coherent rotation mode is degenerate 
with finite q spin waves js^, a circumstance that some- 
times arises due to magnetostatic interactions. The main 
point we wish to raise here is that because of the spin- 
wave Doppler shift, precisely this situation arises when 
the ferromagnet is driven by a current. As in the pre- 
vious subsection, we assume that some type of disorder 
lifts momentum conservation. The effective Hamiltonian 
for the spin waves reads: 

H - ujoblbo + J2 ^i'l)blbq + blJ2 + h-c. (40) 



q=iO 



where is the annihilation operator for the spin wave 
with momentum q and gq- is some unspecified matrix el- 
ement accounting for disorder induced elastic scattering 
of the spin waves. Equation l|40|) is the well Hamilto- 
nian known for a damped harmonic oscillator and can be 
solved exactly. The damping rate for the q = spin wave 
reads: 



^^^^ - T 



dq 



|.gg-p(5(wo - Wg-) 



(41) 



Now we use loq—lo^ = pq^—aq-J^. After a straightforward 
calculation we obtain: 



a\J\ 
47r p2 



(42) 



where we have approximated g^ ~ g. Hence, in the pres- 
ence of elastic spin wave scattering, renormalization of 
the spin wave spectrum due to the current will enhance 



the damping of the lowest spin wave mode. Unlike the 
Gilbert model, the damping rate given by equation H42() 
is independent of wq, implying that the dimensionless 
Gilbert damping coefficient would decline with external 
field if this mechanism were dominant. 



SPIN- WAVE DOPPLER SHIFT AS A 
SPIN-TORQUE EFFECT 

In this section we explain how the influence of an uni- 
form current on magnetization dynamics can be under- 
stood as a special case of a spin-torque eSect^^. The 
latter takes place when a spin current coming from a mag- 
net spin polarized along Mi enters in a magnet spin po- 
larized along M2. In this circumstance there is an imbal- 
ance between the incoming and the outgoing transverse 
component (with respect to M2) of the spin currents in 
magnet 2. Because of spin conservation (resulting from 
the rotational invariance of the system), the imbalance 
in the spin flux across the boundaries of magnet 2 must 
be compensated by a change of the magnetization of that 
magnet, which is described by a new term in the Lan- 
dau Lifshitz equation ^ Q . The microscopic origin of 
the spin current imbalance can be understood as a de- 
structive interference effect, originated by the fact that 
the steady state spin current is a sum over stationary 
states with broad distribution in momentum space |2|. 
Alternatively, it is possible to understand the spin cur- 
rent flux imbalance as a destructive interference in the 
time domain. At a given instant of time, the outgoing 
current-carrying quasiparticles have elapsed a different 
amount of time in magnet 2. This broadening in the in- 
teraction time distribution results in a broadening of the 
spin precession angle [37j |. The average over that distri- 
bution results in a vanishing transverse spin component 
of the outgoing flux. 

The above argument, connecting spin flux imbalance 
and spin-torque, applies to a system in which the inhomo- 
geneous magnetization is described by piecewise constant 
function. It is our contention that the spin wave Doppler 
shift can be understood by applying the same argument 
to the case of smoothly varying magnetization. We con- 
sider again a magnet with charge current j, and spin 
current J7. We assume that current flows in the x di- 
rection and, importantly, that the spin current is locally 
parallel to the magnetization orientation J^{x) = jsQ{x). 
It can be shown that this is the case in a wide range of 
situations. 

The spin density reads S{x) = SoCl{x) where Sq is the 
average spin polarization. We focus on the slab centered 
at X and bounded hy x — dx and x+dx. Spins are injected 
into the slab at the rate js^{x — dx) and leave at the 
rate js^{x -\- dx). The resulting spin current imbalance 
is 2dxjsdx^. Therefore, there must be a spin transfer to 
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the local magnetization: 
dS{x) 



dt 



(43) 



ST 



Now, using \ (}\^ = 1 at every point of the space we obtain: 



dSix) 



dt 



■ js^{x) X (d^nix) X n{x)) (44) 



ST 



which is exactly the same result obtained in|T3. Includ- 
ing this term in the Landau Lifshitz equation and solving 
for small perturbations around the homogeneous ground 
state (spin waves) results into the spin wave Doppler shift 
discussed in previous sections. In conclusion, this argu- 
ment demonstrates that the spin- wave Doppler shift and 
spin transfer torques are different limits of the same phys- 
ical phenomena, the transfer of angular momentum from 
the quasiparticles to the collective magnetization when- 
ever the latter is not spatially uniform. 

DISCUSSION AND CONCLUSIONS 

The effect of high current densities on the magneti- 
zation dynamics of ferromagnetic metals have been ex- 
plored experimentally in several configurations. In point 
contact experiments, a large current density is injected 
from a normal metallic contact into a ferromagnetic mul- 
tilayer 0, or single layer When a large flow of 
electrons ( current density j~ 10* A cm~^) enters into 
the ferromagnetic multilayer, the resistivity presents an 
abrupt increase which has been related to the coherent 
precession of spin waves I3| and/or phonons The 
fact that Ji and Chien [8| report similar results when 
the current is injected into a a single ferromagnetic layer 
demonstrates that interlayer coupling is not essential for 
the anomalies observed in transport. It must be noted 
that when the current flow is such that the electrons go 
from the ferromagnetic layer(s) toward the point contact, 
no anomaly is observed. Similar transport anomalies at 
currents densities higher than those of current induced 
magnetization switching are observed by a number of dif- 
ferent groups in a system of two adjacent ferromag- 
netic nanopillars. In this system a large current density 
flows from one ferromagnet to the other. 

The fact that the current densities at which the anoma- 
lous behavior takes place is of the same order of magni- 
tude than the current at which the spin wave Doppler 
shift makes the collinear state unstable might lead to 
suggest a connection between the two. However, the ex- 
periments in the point contact geometry show that the 
transport anomalies only occur for one direction of the 
current, something which seems at odds with the spin 
wave Doppler shift instability. 



In summary, the focus of this paper is on the effect of 
the current in the spin wave dynamics of a bulk ferro- 
magnetic metal. We have addressed two types of effects: 
the change in the spin wave dispersion and change in 
the spin wave damping. These quantities are given, at 
a formal level, by the spin wave propagator. The cen- 
tral idea is that the spin wave propagator is a functional 
of the quasiparticle occupation function. In the pres- 
ence of the current the occupation function changes, af- 
fecting both the dispersion and the damping of the spin 
waves. Throughout the paper we have assumed that the 
functional relation between the quasiparticle occupation 
function and the spin wave propagator remains the same 
when the system is out of equilibrium. In that sense, 
the above derivations are heuristic. Our main conclu- 
sions are: i) A current j flowing through a metallic ferro- 
magnet results in a spin current S which modifies its 
spin wave spectrum by an amount proportional to q-J^. 
a) This modification, which was derived by BJZ for a 
fully polarized s-d model, occurs as well in a non fully 
polarized Hubbard model, in which the d electrons are 
itinerant and, according to the arguments of section II, 
in typical real- world ferromagnets. in) In the presence 
of elastic two magnon scattering, the spin wave Doppler 
shift leads as well to a broadening of the lowest spin wave 
mode (Eq. I42II . which is proportional to the spin current. 
iv) Both the spin- wave Doppler shift in spatially homo- 
geneous ferromagnets and the spin torque effect in inho- 
mogenoeus structures 0, are a consequence of the spin 
transfer from the quasiparticles to the collective magne- 
tization when the latter is spatially inhomogeneous. 

We acknowledge fruitful discussions with M. Abolfath 
and M. Tsoi. Work at the University of Texas was sup- 
ported by the Welch Foundation and by the National Sci- 
ence Foundation under grants DMR-0210383 and DMR- 
0115947. Work at University of Alicante supported by 
Ramon y Cajal Program, Ministerio de Ciencia y Tec- 
nologia. 



[2] 
[3] 
[4] 

[5] 
[6] 
[7] 



Spin Dependent Transport in Magnetic Nano structures, 
Ed. by S. Maekawa and T. Shinjo, (Taylor and Francis, 
2002). 

J.C. Slonczewski, J. Mag. Mat. Mag. 159, LI (1996). 
L. Berger, Phys. Rev. B54, 9353 (1996). 
M. Tsoi et al, Phys. Rev. Lett. 80, 4281 (1998); M. Tsoi 
at al.. Nature 406, 46 (2000); 

M. Tsoi et al, Phys. Rev. Lett. 89, 246803 (2002). 
J.Z. Sun, J. Magn. Mag. Mater. 202, 157 (1999). 
E. B. Myers et al. Science 285, 867 (1999) J.A. Katine et 
al, Phys. Rev. Lett. 84, 4212 (2000); E.B. Myers, et al., 
Phys. Rev. Lett. 89, 196801, (2002); S.I. Kiselev em et 
al., preprint cond-mat/0306259 (2003). W.H. Rippard, 
M.R. Pufall,and T.J. Silva, Appl. Phys. Lett. 82, 1260- 
1262 (2003). F. B. Mancoff , et al Appl. Phys. Lett. 83, 
1596 (2003). 



12 



[8] Y. Ji and C.L. Chien and M. D. Stiles, Phys. Rev. Lett. 

90, 106601 (2003). 
[9] S. Urazhdin et al., Phys. Rev. Lett. 91, 146803 (2003) 

[10] J. Z. Sun, Phys. Rev. B 62, 570 (2000 ); A. Brataas, 
Y. V. Nazarov, and G.E.W. Bauer, Phys. Rev. Lett. 84, 
2481 (2000); X. Waintal and P.W. Brouwer, Phys. Rev. 
B 63, 220407 (2001); C. Heide, Phys. Rev. B 65, 054401 
(2002); M. Stiles and A. ZangwiU, Phys. Rev. B 65, 
014407 (2002); J.-E. Wegrowe, Appl. Phys. Lett. 80, 3775 
(2002); S. Zhang, P.M. Levy, and A. Pert, Phys. Rev. 
Lett. 88, 236601 (2002); G.E.W. Bauer, Y. Tserkovnyak, 
D. Huertas-Hernando, and A. Brataas, Phys. Rev. B 
67, 094421, (2003); M.L. Pohanski and P.W. Brouwer, 
preprint cond-mat/0304069 (2003); A. Shapiro, P. M. 
Levy, and S. Zhang Phys. Rev. B 67, 104430 (2003); A. 
Pert et al, preprint cond-mat/0310737 (2003). 

[11] Y.B. Bazaliy, B. A. Jones and S.-C. Zhang, Phys. Rev. 
B57, R3213 (1998) 

[12] More precisely, J = ^J^k ^ [^l ' ^i]' where TV is 
the number of sites in the lattice, and | and J. are defined 
in the axis of the average magnetization. 

[13] C. Herring and C. Kittel, Phys. Rev. 81, 869 (1951). 

[14] E.P. Wolfarth, Rev. Mod. Phys. 25, 211 (1953); D.A. 
Papaconstantopoulos, Handbook of the Band Structure 
of Elemental Solids (Plenum Press, New York, 1986). 

[15] R. E. Prange and V. Korenman, Phys. Rev. B 19, 4691 
(1979) 

[16] T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 
netism, (Springer- Verlag, Berlin, 1985) 
[17] J. Hubbard, Phys. Rev. B19, 2626 (1979) 
[18] H.J. Schulz, Phys. Rev. Lett65, 2462 (1990) 
[19] Sections 3.3 and 3.4 in E. Fradkin, Field Theories of Con- 
densed Matter Systems, (Addison Wesley, Reading, MA, 
1991) 



[20] S. Doniach and E. H. Sondheimer, Green's Functions for 
Solid State Physicists, W. A. Benjamin, Inc (1974) 

[21] R.L. Stratonovitch, Dok. Akad. Nauk. SSSR 115, 1097 
(1957); [Sov. Phys. Dokl. 2, 416 (1958)]. J. Hubbard, 
Phys. Rev. Lett. 3, 77 (1959) 

[22] J.W. Negele and H. Orland, Quantum Many-Parttcle 
Systems, Addison Wesley, (1988) 

[23] Assa Auerbach, Interacting Electrons and Quantum Mag- 
netism, Springer- Verlag, 1994. 

[24] R. Skomski and J. M. D. Coey, Permanent Magnetism 
(Institute of Physics, Bristol, 1999). 

[25] Claude Cohen- Tannoudji, Gilbert Grynberg, Jacques 
Dupont-Roc Atom-Photon Interactions: Basic Processes 
and Applications, (Wiley-Interscience, New York, 1998) 

[26] V. Korenman, J. L. Murray, and R. E. Prange Phys. Rev. 
B 16, 4032 (1977) 

[27] A. J. Millis, P.B. Littlewood and B. I. Shraiman, Phys. 
Rev. LettT4, 5144 (1995) 

[28] P. M. Levy and S. Zhang, Phys. Rev. Lett. 79 , 5110 
(1997) 

[29] G. Tatara and H. Fukuyama, Phys. Rev. Lett. 78, 3773 
(1997) 

[30] A. Pert, J. Phys. C. 2, 1784 (1969) and references therein. 
[31] V. Kambersky, Canadian Journal of Physics, 48 2906 
(1970) 

[32] V. Kambersky and C. E. Patton, Phys. Rev. B 11, 2668 
(1975) 

[33] V. Korenman, J. L. Murray, and R. E. Prange Phys. Rev. 

B 16, 4048-4057 (1977) 
[34] L. Berger, J. Phys. Chem. Solids 38, 1321 (1977) 
[35] N. W. Ashcroft and N. D. Mermin, Solid State Physics, 

W.B. Saunders Company (1976) 
[36] R. Arias and D. Mills, Phys. Rev. B 60, 7395 (1999) 
[37] J. Fernandez-Rossier et al, cond-mat/0304492 (2003) 



